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ABSTRACT 

Context. NGC 288 is a globular cluster with a well developed blue horizontal branch covering the so-called w-jump which indicates the onset 
of diffusion. It is therefore well suited to study the effects of diffusion in blue horizontal branch (HB) stars. 

Aims. We compare observed abundances to predictions from stellar evolution models calculated with diffusion and from stratified atmospheric 
models. We verify the effect of using stratified model spectra to derive atmospheric parameters. In addition we investigate the nature of the 
overluminous blue HB stars around the zz-jump. 

Methods. We define a new photometric index sz from uvby measurements that is gravity sensitive between 8 000K and 12 000 K. Using 
medium-resolution spectra and Stromgren photometry we determine atmospheric parameters (J e g, log g) and abundances for the blue HB stars. 
We use both homogeneous and stratified model spectra for our spectroscopic analyses. 

Results. The atmospheric parameters and masses of the hot HB stars in NGC 288 show a behaviour seen also in other clusters for temperatures 
between 9 000K and 14 000 K. Outside this temperature range, however, they follow rather the results found for such stars in tuCen. The 
abundances derived from our observations are for most elements (except He and P) within the abundance range expected from evolutionary 
models that include the effects of atomic diffusion and assume a surface mixed mass of 10~ 7 M 0 . The abundances predicted by stratified model 
atmospheres are generally significantly more extreme than observed, except for Mg. The use of stratified model spectra to determine effective 
temperatures, surface gravities and masses moves the hotter stars to a closer agreement with canonical evolutionary predictions. 

Conclusions. Our results show definite promise towards solving the long-standing issue of surface gravity and mass discrepancies for hot HB 
stars, but there is still much work needed to arrive at a self-consistent solution. 
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1. Introduction 

Low-mass stars burning helium in a core of about 0.5 M G and 
hydrogen in a shell populate a roughly horizontal region in the 
optical colour-magnitude diagrams of globular clusters, which 
has earned them the name “horizontal-branch” (HB) stars (ten 
Bruggencate, 1927). The hot (or blue) HB stars near an effec- 
tive temperature of 1 1 500 K are of special interest because they 


* Based on observations with the ESO Very Large Telescope at 
Paranal Observatory, Chile (proposal ID 71.D-0131) 


exhibit a number of intriguing phenomena associated with the 
onset of diffusion. 


1.1. Diffusion 

A large photometric survey of many globular clusters by 
Grundahl et al. (1999) demonstrated that the Stromgren u- 
brightness of blue HB stars suddenly increases near 11 500 K. 
This so-called w-jump is attributed to a sudden increase in the 
atmospheric metallicity of the blue HB stars to super-solar val- 
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ues, due to the radiative levitation of heavy elements. A simi- 
lar effect can be seen in broad-band U,U — V photometric data 
(Ferraro et al. 1998, Gl). The work of Behr et al. (1999, 2000b) 
and Moehler et al. (2000) confirmed with direct spectroscopic 
evidence that the atmospheric metallicity does indeed increase 
to solar or super-solar values for HB stars hotter than the u- 
jump. A list of earlier observations of this effect can be found 
in Moehler (2001, Sect. 2.4). Later studies include Behr (2003, 
M3, M 13, M 15, M 68, M92, and NGC 288), Fabbian et al. 
(2005, NGC 1904), and Pace et al. (2006, NGC 2808). These 
findings also helped to understand the cause of the so-called 
“low-gravity” problem: Crocker et al. (1988) and Moehler et 
al. (1995, 1997) found that hot horizontal branch stars - when 
analysed with model spectra of the same metallicity as their 
parent globular cluster - show significantly lower surface grav- 
ities than expected from evolutionary tracks. Analysing them 
instead with more appropriate metal-rich model spectra re- 
duces the discrepancies considerably (Moehler et al. 2000). 
The more realistic stratified model atmospheres of Hui-Bon- 
Hoa et al. (2000) and LeBlanc et al. (2009) reduce the discrep- 
ancies in surface gravity even further (see below for more de- 
tails). Along with the enhancement of heavy metals a decrease 
in the helium abundance by mass Y is observed for stars hotter 
than sail 500 K, while cooler stars have normal helium abun- 
dances within the observational errors. The helium abundance 
for these hotter stars is typically between 1 and 2 dex smaller 
that the solar value (e.g. Behr 2003). A trend of the helium 
abundance relative to T e ff was discussed by Moni Bidin et al. 
(2009, 2012). Finally, blue HB stars near »11 500 K show a 
sudden drop in their rotation rates (Peterson et al. 1995, Behr 
et al. 2000a and Behr 2003), and in some globular clusters (e.g., 
M 13) a gap in their HB distribution. 

The possibility of radiative levitation of heavy elements and 
gravitational settling of helium in HB stars had been predicted 
long ago by Michaud et al. (1983), but the discovery of its very 
sudden onset near 1 1 500 K was a complete surprise. Quievy et 
al. (2009) have shown that helium settling in HB stars cooler 
than ~1 1 500 K is hampered by meridional flow. In stars hot- 
ter than this threshold, helium may settle and the superficial 
convection zone disappears and therefore diffusion may come 
at play in superficial regions of these stars. Recent evolution- 
ary models of HB stars that include atomic diffusion calculated 
by Michaud et al. (2008, 2011) can reproduce the abundance 
anomalies of several elements observed in blue HB stars rea- 
sonably well. However, these models do not treat the atmo- 
spheric region in detail. Instead they assume an outer superfi- 
cial mixed zone of approximately 1O~ 7 M 0 below which sepa- 
ration occurs. 

The detection of vertical stratification of some elements, 
especially iron, in the atmospheres of blue HB stars lends 
additional support to the presence of atomic diffusion there 
(Khalack et al. 2007, 2008, 2010), but is at variance with the 
mixed zone introduced by Michaud et al. (2008, 2011). Hui- 
Bon-Hoa et al. (2000) and LeBlanc et al. (2009) present stellar 
atmosphere models of blue HB stars which include the effect of 
vertical abundance stratification on the atmospheric structure. 
These models estimate the vertical stratification of the elements 
caused by diffusion by assuming that an equilibrium solution 


(i.e. giving a nil diffusion velocity) for each species is reached. 
Assuming a sudden onset of atomic diffusion near 11 500K, 
these models predict photometric jumps and gaps (Grundahl et 
al. 1999, Gl in Ferraro etal. 1998) consistent with observations 
(see LeBlanc et al. 2010 for more details). The photometric 
changes, as compared to chemically homogeneous atmosphere 
models, are due to the modification of the atmospheric struc- 
ture caused by the abundance stratification. 

1.2. NGC 288 

The z/-jump described above can be clearly seen in the colour- 
magnitude diagram of NGC 288 shown in Fig. 1, where one 
also finds a group of stars with large (and unexplained) scatter 
in their u magnitudes and (u-y) colours around the jump region 
(triangles in Fig. 1). With maximum errors of 0 1 1'008 in u and 
0 1 ? 1 003 in y it is unlikely that their positions are caused by pho- 
tometric errors. The evolutionary status of these stars is unclear. 
While their bright u magnitudes are suggestive of radiative lev- 
itation, the effect would have to be extreme and some of them 
appear to lie on the cool side of the w-jump. Similar groups of 
unexpectedly bright stars can be found in other globular clus- 
ters (see M 2 and M 92 in Grundahl et al. 1999), with NGC 288 
presenting the most pronounced case. We refer the reader to 
Sect. 8 for a more detailed discussion. 

A colour spread along the red giant branch in NGC 288 first 
reported by Yong et al. (2008) was identified as a split by Roh 
et al. (20 11), which was confirmed by Carretta et al. (20 1 1 ) and 
Monelli et al. (2013). In their excellent review on second and 
third parameters to explain the horizontal-branch morphology, 
Gratton et al. (2010) estimated that a range in helium abun- 
dance of AT = 0.012 would explain the temperature range of 
the horizontal branch in NGC 288. This is consistent with the 
helium range found more recently from the analysis of main 
sequence photometry by Piotto et al. (2013). Such a small vari- 
ation in helium is unfortunately too small to be detected by our 
analysis. 

2. Observations 

The targets were selected from the Stromgren photometry of 
Grundahl et al. (1999). We selected 71 blue horizontal-branch 
star candidates (see Fig. 1) and 17 red giants. Of the blue HB 
candidates three turned out to be red HB stars and one had ex- 
tremely noisy spectra. In this paper we discuss only the ob- 
servations of the 67 bona-fide hot horizontal-branch stars (see 
Table 1 for their coordinates and photometric measurements). 

The spectroscopic data were obtained between July 3 and 
27, 2003 (date at the beginning of the night, see Table 2 for 
details) in Service Mode using the multi-object fibre spectro- 
graph FLAMES+G1RAFFE (Pasquini et al. 2000), which is 
mounted at the UT2 Telescope of the VLT. The fibre systems 
MEDUSA1 and MEDUSA2 allow to observe up to 132 ob- 
jects simultaneously. We used the low spectroscopic resolution 
mode with the spectral ranges 3620 - 4081 A (LR1, T/AT = 
8000), 3964 - 4567 A (LR2, T/A A = 6400) and 4501 - 5078 A 
(LR3, T/A/l = 7500). GIRAFFE had a 2kx4k EEV CCD chip 
( 15/um pixel size), with a gain of 0.54 e~ ADLL 1 and a read-out- 
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Fig. 1. The blue HB stars in NGC 288 as observed by Grundahl et al. (1999). The subsample, for which GIRAFFE spectra were observed, is 
marked by filled circles. The dot-dashed line in the left part marks the u-y colour of the w-jump. Stars on the red side of the line do not show 
evidence for radiative levitation, while stars on the blue side do. The triangles mark overluminous stars around the (/-jump (from red to blue: 
122, 103, 127, 101, 146, 142, 100, 183). 


noise of 3.2 e - . Each night four screen flat fields, five bias, and 
one ThAr wavelength calibration frame were observed. In ad- 
dition to the daytime ThAr spectra covering all fibres, we also 
observed simultaneous ThAr spectra during the science obser- 
vations in five of the fibres. Dark exposures were not necessary, 
because the dark current of 2 e - h 1 pixel -1 is negligible. 

3. Data reduction 

The spectroscopic data were reduced using the girBLDRS soft- 
ware (http://girbldrs.sourceforge.net/, version 1.10) and ESO 
MIDAS (see Drews 2005 for details). The 2-dimensional bias 
and flat field frames were averaged for each night. The aver- 
aged flat fields were used to determine the positions and widths 
of the fibre spectra. 64-pixel portions of each fibre spectrum 
were averaged and fitted with a point spread function (PSF). 
A polynomial fit to the PSF parameters then provides the PSF 
for the whole frame, which is used for all extractions later on. 
The one-dimensional flat field spectra extracted using Home’s 


method (Horne 1986) are the so-called “narrow flat-fields” 
(nff), which are used to correct the extracted science spectra 
for spectrograph and detector signature. As the nff spectra are 
not normalised they still show the spectral signature of the flat 
field lamp, which however varies slowly with wavelength. 

The daytime Th-Ar wavelength calibration spectra were ex- 
tracted for each fibre via a simple sum, the spectral lines lo- 
calised and fit by an analytical model. The presence of simulta- 
neous ThAr spectra in the science frames allowed us to refine 
the localisation and wavelength calibration of the science data 
by correcting for residual wavelength shifts between daytime 
and nighttime observations. The refinement was repeated un- 
til the difference between the observed ThAr positions and the 
laboratory ones was below 0.001 kms -1 . The wavelength off- 
sets obtained from the cross-correlation are then linearly inter- 
polated across the CCD (spatial direction), and the interpolated 
offsets are applied to each fibre. 

The raw science spectra were bias and flatfield corrected 
and finally rebinned to constant wavelength steps. The median 
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of the sky signal obtained in 1 8 dedicated fibres was subtracted 
from the science data. 

The spectra for a given star and setup were then normalised 
as described in Moehler et al. (2011) taking care to use only re- 
gions free from strong lines for the continuum definition. This 
normalisation allows us to exclude outlying pixels during the 
averaging of the spectra. 


4. Radial velocities 

After the barycentric correction, the observed spectra were first 
co-added (without further radial velocity correction) and fitted 
with stellar model spectra to get a first estimate of their effective 
temperatures, surface gravities, and helium abundances (see 
Sect. 5.2 for details). The individual spectra of each star were 
then cross-correlated (see Tonry & Davis 1979 for more de- 
tails) with the best fitting synthetic spectrum derived this way. 
Only regions of hydrogen or helium lines were selected prior 
to the cross-correlation. The peak of the cross-correlation func- 
tion was fitted with a Gaussian function to determine the ra- 
dial velocity to sub-resolution accuracy. The velocity-corrected 
spectra were co-added and fitted with synthetic model atmo- 
spheres (see Sect. 5). In a second step, the best-fit synthetic 
spectra were then used to repeat the cross-correlation. The 
lcr errors in radial velocity for individual spectra range from 
1.8kms _1 to up to lOkms -1 depending on the quality. The me- 
dian error is 2.7 km s~ ! . 

For each star, we calculated the standard deviation of the ra- 
dial velocity measurements of the individual spectra. We then 
counted, how many of the individual measurements deviate 
more than lcr from the mean. In none of the stars this was larger 
than 30%, i.e. the scatter is consistent with a random distribu- 
tion within the measurement error. We thus find no evidence 
for radial velocity variations that might indicate the presence 
of close binaries in our sample. We note that this absence of 
evidence should not be confused with an evidence of ab- 
sence, i.e. we cannot say anything about the frequency of 
binaries among our target stars, as neither our method nor 
our observations were tailored towards the detection of bi- 
naries. There is also no evidence for field contamination in our 
sample. 

The radial velocities of the 67 stars range from 
-48.6 km s~' to -36.9 km s~' with a median value of 
-42.6+2.7 km s _1 , that is close to the average value of 
-45.4 km s -1 given by Harris (1996) and to the recently pub- 
lished values of -43.5 km s~' (Szekely et al. 2007, uncertainty 
1-2 km s^ 1 ), -46. 15+2.55 km s~' (Carretta et al. 2009a), and 
-45.1+0.2 km s _1 (Lane et al. 2010). 

Each spectrum was corrected to laboratory wavelength. As 
an additional safeguard against outliers only a subset of pixel 
values around the median (11 out of 13 for LR1, 5 out of 7 
for LR2, 3 out of 5 for LR3) was averaged. In order to allow a 
combination of the three wavelength ranges into one spectrum 
for fitting the LR1 data were convolved with a Gaussian with 
a FWHM of 0.48 A to have the same resolution for all three 
regions. 


5. Atmospheric Parameters 

5.1. Determination from photometric data 

We used the Stromgren photometry of Grundahl et al. 
(1999) to get a first estimate of the effective tempera- 
tures and surface gravities of the stars redwards of the u- 
jump, similar to our work in Moehler et al. (2003). As ref- 
erence we used theoretical colours from Kurucz (ATLAS9, 
1993) for metallicities [M/H] = -1.0 and -1.5. The metal- 
licity [Fe/H] of NGC 288 given by Carretta et al. (2009b) is 
-1.32. Unfortunately we have no photometry for the stars. 
Therefore, guided by the equations of Moon & Dworetsky 
(1985) we looked for a combination of colours that provides 
a rectangular grid in the 7% , log g plane for the range 8 000 K 
L T e ff < 12 000 K. We used the definition for a 0 from Moon & 
Dworetsky, but we defined in addition an index sz as follows 

«o = 1-36 • ( b —y) + 0.36 ■ m\ + 0.18 • c\ - 0.2448 

sz = —0.07 • (b - y) — m\ + 0.1 • c\ + 0.1 

As one can see from Fig. 2 ao correlates with 7% and sz 
correlates with log g and T e ff. The photometric data were dered- 
dened with 7 ?b-v = O'.' 1 03 (Carretta et al. 2000), using 7?b- y = 
0.75 -£ B -v- 



Fig. 2. Theoretical ao, sz grid as derived from Kurucz (1993) colours 
for [Fe/H] = -1.5. In addition we show the positions of the cool blue 
HB stars (i.e. redward of the vertical line in Fig. 1) in this diagram as 
observed by Grundahl et al. (1999). 

To restrict the range of fitting we used only logg between 
3 and 4 and T eS between 8 000 K and 12 000 K. We fitted a 2 nd 
order polynomial to the relation 7% (ao), which yielded an rms 
deviation of 30 K. For the surface gravities we fitted 2 nd order 
polynomials to the relation logg(ao, sz), which yielded an rms 
deviation of 0.04 in logg. Temperatures and surface gravities 
derived from these relations are listed in Table 3 and plotted in 
Fig. 3 for stars cooler than the z/-jump. 

The errors provided in Table 3 are a quadratic combination 
of errors from the photometric data, the uncertainty of the red- 
dening in the y band (estimated to be 0™02) and the error of the 
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fit to the theoretical relations. In order to have an estimate of the 
influence of metallicity we compared the results obtained for 
the two metallicities mentioned above. The differences in tem- 
perature vary almost linearly from -0.1% at 8000 K to +0.8% 
at 12 000 K, with the results from the more metal-poor models 
being hotter at the hot end. The differences are therefore well 
below the errors on the individual temperatures. 

The differences in surface gravities for these two metallici- 
ties show a quadratic dependency on the effective temperature, 
varying from +0.02 at about 8 600 K (with the surface gravities 
from the metal-poor models being higher at 8 600 K than those 
from the more metal-rich models) to 0 at 9 000 K and 1 1 000 K 
via -0.02 at about 10 000 K. For temperatures below 8 500K 
the differences increase to -0.08, but those stars are too cool to 
be analysed spectroscopically with our methods (see Sect. 5.2). 
Within the temperature range of interest for further analysis the 
differences are below the average errors of the surface gravities 
given in Table 3. For reasons of consistency between photomet- 
ric and spectroscopic analysis we decided against averaging the 
photometric grids. 

Talcing together the errors listed in Table 3 and the uncer- 
tainties due to the metallicity, we estimate typical errors in Teff 
and logg to be about 5% and 0.13 dex, respectively. 



Fig. 3. Atmospheric parameters derived from the a 0 ,sz indices for 
stars not showing evidence for radiative levitation. For comparison we 
also show curves representing a canonical zero-age horizontal branch 
and terminal age horizontal branch from Moehler et al. (2003). 


5.2. Line profile fitting with homogeneous model 
spectra 

For the stars redder than the w-jump, we computed model at- 
mospheres with convection for [M/H] = -1.5 using ATLAS9 
(Kurucz 1993) and used Lemke’s version 1 of the LINFOR 
program (developed originally by Holweger, Steffen, and 

1 For a description see http://a400.stemwarte.uni- 

erlangen.de/~ai26/linfit/linfor.html 


Steenbock at Kiel University) to compute a grid of theoreti- 
cal spectra which include the Balmer lines H,, to H 22 , Fie 1 and 
He 11 lines. The grid covers the range 7 500 K < r e ff < 12 000 K, 
2.5 < logg < 5.0, and -3.0 < log ^ < -1.0. For the stars 
bluer than the w-jump we used model atmospheres computed 
for [M/H] = +0.5 (see Moehler et al. 2000 for details). From 
an extrapolation of the LTE/NLTE threshold for subdwarf B 
stars (Napiwotzki 1997) we assume that LTE is a valid approx- 
imation here as well. To establish the best fit to the observed 
spectra, we used the fitsb2 routines developed by Napiwotzki 
et al. (2004), which employ aqf 2 test. The cr necessary for the 
calculation ofx 2 is estimated from the noise in the continuum 
regions of the spectra. The fit program normalises synthetic 
model spectra and observed spectra using the same points for 
the continuum definition. 

While the errors were obtained via boot strapping and 
should therefore be rather realistic, they do not include possible 
systematic errors due to, e.g., flat field inaccuracies or imper- 
fect sky subtraction. The true errors in T e ff are probably close 
to those from photometry, i.e. 5%, and the true errors in logg 
are probably about 0.1. 

Because the x 2 fit °f the line profile can lead to ambigu- 
ous results for effective temperatures close to the Balmer max- 
imum, the results from Sect. 5.1 are taken as initial parameters 
for the spectral line profile fitting procedure. 

In Table 3 we list the results obtained from fitting the 
Balmer lines H /; to H 10 (excluding H, to avoid the Can H line) 
in the cool stars. For the cool stars we did not fit the helium 
lines, as they are very weak and tend to produce spurious re- 
sults. We verified, however, that the helium lines predicted by 
the model spectra with solar helium abundance did reproduce 
the observations. The reddest stars showing many strong metal 
lines (55, 60, 61, 63, 72) were omitted from the line profile fits 
because our model spectra contain only H and He lines. For the 
stars hotter than 1 1 500 K we fitted in addition the He 1 lines at 
AA 4026A, 43 8 8 A, 4472A, and 4922A. 

Fig. 4 shows a gap at 8 500K< T e s <9 000K in the at- 
mospheric parameters derived from line profile fitting, that was 
first noted and discussed by Moehler et al. (2003) for obser- 
vations of HB stars in M 13. The stars populating this region 
in Fig. 3 close to the zero-age horizontal branch (ZAHB) are 
in Fig. 4 at about 8 000 K, above the terminal-age horizontal 
branch (TAHB). Prompted by the referee we treated a homo- 
geneous model spectrum calculated for T eS - 8 800 K, logg = 
3.35 and [M/H] = -1.5 as described in Sect. 6.3 to simulate an 
observed spectrum. Varying the starting values of the fit from 
8 000 K to 9 500 K for T e t f and 3.2 to 3.6 for logg did not affect 
the fit results which stayed at T e s — 8 841 K and logg = 3.35. 
Varying the FWHM of the Gaussian, with which the model 
spectra are convolved before fitting, from 0.6 A to 0.75 A re- 
sulted in T efs between 8 818 K and 8 862 K and logg between 
3.34 and 3.36. We noted, however, that the resolution of the 
observed data seems to vary slightly with wavelength — at the 
blue end the predicted line profiles are very slightly wider and 
deeper than the observed ones, while for ¥tp the observed line 
core is slightly more narrow and deep. To rule out possible 
small variations in resolution across the wavelength range as 
the cause for the gap, we convolved the observed spectra to 
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Table 4. Temperatures, surface gravities, and helium abundances from 
line profile fits for stars bluer than the (/-jump and showing evidence 
of radiative levitation). The errors are derived from the y 2 fit (see 
Sect. 5.2 for details) and are statistical errors only. In addition we list 
the helium abundance as determined with LINFOR (see Sect. 6 for 
details). 


ID 

Teff 

logg 

log(— ) log(— ) 

b V «H / fitsb2 ® V «H M.TNFOR 

100 

1 1400±100 

3.65±0.03 

— 2.23±0. 1 1 

-2.14 

101 

1 1400±100 

3.78±0.03 

— 1 .95±0. 1 1 

-1.93 

142 

1 1 400± 1 00 

3.72±0.03 

-2.98±0.12 

-2.76 

146 

1 1 600± 1 00 

3.81±0.03 

-2.06±0.10 

-2.02 

179 

1 1 600± 1 00 

3.77±0.03 

-2.37±0.10 

-2.14 

183 

12300±100 

3.93±0.03 

-2.27±0.12 

-2.36 

187 

1 1 700± 100 

3.84±0.03 

-2.74±0.13 

-2.87 

195 

12400 ±100 

3.86±0.03 

-2.55±0.15 

-2.51 

196 

12300±100 

3.91±0.03 

-2.67±0.18 

-2.49 

199 

12400 ±100 

3.94±0.03 

-2.79±0.18 

-2.84 

212 

12900 ±100 

3.93±0.02 

-2.66±0.12 

-2.57 

213 

131 00± 1 00 

3.97±0.03 

— 2.61 ±0. 1 5 

-2.68 

216 

12900 ±100 

3.94±0.02 

-2.66±0.14 

-2.71 

221 

13200±100 

3.97±0.03 

-2.67±0.13 

-2.60 

228 

13800±100 

4.00±0.03 

-2.43±0.17 

-2.43 

230 

13800±100 

4.06±0.03 

-2.28±0.10 

-2.34 

231 

13300±100 

4.00±0.02 

-2.72±0.09 

-2.55 

240 

13300±100 

3.99±0.03 

— 2.7 1 ±0. 1 1 

-2.61 

242 

14300±100 

4.03±0.03 

-2.37±0.13 

-2.28 

243 

14400 ±100 

4.09±0.04 

-2.52±0.15 

-2.43 

275 

15400±200 

4.16±0.03 

-2.44±0.15 

-2.27 

288 

15200±200 

4.14±0.04 

— 2. 19±0. 1 1 

-2.16 

292 

151 00± 1 00 

4.23 ±0.03 

-2.92±0.07 

-2.80 

300 

15400±200 

4.14±0.04 

-2.21±0.10 

-2.19 

304 

16900±200 

4.29±0.04 

-2.57±0.12 

-2.46 

347 

16400±200 

4.38±0.04 

-2.11 ±0.09 

-2.06 


a resolution of 2.5 A, which is close to the resolution of the 
data of Moni Bidin et al. (2007, 2009, 2011) and Salgado et al. 
(2013), which do not show such a gap in their results although 
they used the same analysis methods as we do. Unfortunately 
the gap remains in the results of the line profile fitting. Also 
using model atmospheres with or without convection did not 
affect the gap. 

The stars in Fig. 4 can be split in two groups with respect 
to their position relative to the evolutionary sequences: The 
stars with effective temperatures between 9 000 K and roughly 
14 000 K lie between the ZAHB and the TAHB, while the 
stars hotter or cooler than this temperature range lie above the 
TAHB. 

In order to see how the parameters derived for the hot HB 
stars in NGC 288 compare to those obtained for other clus- 
ters and to the theoretical predictions we calculated the dif- 
ference between the derived surface gravity and the one on 
the ZAHB for the given temperature. Fig. 5 compares our re- 
sults to those obtained for M 80, NGC 5986 (Moni Bidin et al. 
2009), NGC 6752 (Moni Bidin et al. 2007), M22 (Salgado et 
al. 2013), and wCen (Moehler et al. 2011; Moni Bidin et al. 
2011). Here the NGC 288 results look like those from wCen 
for temperatures below 9 000 K or above about 14 000 K, while 



Fig. 4. Atmospheric parameters derived from the line profile fitting for 
the stars hotter (full symbols) and cooler (open symbols) than the u- 
jump. For comparison we also show a canonical zero-age horizontal 
branch and temiinal age horizontal branch from Moehler et al. (2003). 


they coincide with the other clusters for temperatures between 
9 000 K and 14 000 K. 

Using the spectroscopically determined atmospheric pa- 
rameters, y magnitudes and the distance modulus of the globu- 
lar cluster (( m - A7) v = 14’"95, Carretta et al. 2000), we derived 
the masses shown in Fig. 6 using the equation 

log A i_ = 36 ; 10 7 . ]Q-0My-(m-M)r-y t k) 

M Q n ■ g G • R 0 [pc] 

where y^ is the theoretical v magnitude at the stellar surface 
from Kurucz (1993) and R 0 is the solar radius in parsec. 
Obviously the masses are systematically too low, except for 
the stars just hotter than the z/-jump. When comparing to re- 
sults from other clusters we see that the masses of the stars in 
NGC 288 show a similar behaviour as the surface gravities with 
respect to HB stars in other globular clusters. 

6. Abundances 

We used the parameters derived in Sect. 5.2 for the follow- 
ing abundance analysis. 

6.1. Helium 

Helium abundances were already derived during the determi- 
nation of effective temperatures and surface gravities. In addi- 
tion, we determined them (together with abundances for other 
elements) via spectrum synthesis using the abundance fitting 
routine of LINFOR (see Sect. 6.2 for details). The resulting 
abundances are listed in Tables 3 and 4. For the cool stars the 
abundances should be treated with caution as the helium lines 
are rather weak. However, the average fitted helium abundance 
log ^ of -1.03 agrees well with our assumption of solar he- 
lium abundance for the cool stars. For the hot stars the helium 
abundance derived with LINFOR are in general slightly higher 
than the ones derived with fitsb2, but in 80% of the cases the 
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Fig. 5. Surface gravity relative to the zero-age horizontal branch at the 
given effective temperature. For comparison we also show a canoni- 
cal zero-age and terminal-age horizontal branch from Moehler et al. 
(2003). In the upper plot we show also results from FORS2 obser- 
vations of hot HB stars in M 80, NGC 5986 (Moni Bidin et al. 2009), 
NGC 6752 (Moni Bidin et al. 2007), and M 22 (Salgado et al. 2013). In 
the lower plot we also provide the results obtained from FLAMES and 
FORS2 observations of hot horizontal branch stars in wCen (Moehler 
et al. 2011, small squares; Moni Bidin et al. 2011, small triangles). 


difference is below the error provided by fitsb2. On average the 
LINFOR abundances are higher by 0.06 dex, i.e. well below the 
error of the fitsb2 results. 

As helium abundances have been derived for HB stars in 
many globular clusters with similar methods and model spec- 
tra as used in this paper, we compare our results to published 
data in Fig. 7. The helium abundances of Moni Bidin et al. 
(2007, 2009, 2012) and Moehler et al. (2000, 2003) were de- 
rived from low resolution data with 2.6 A- 3.4 A resolution. 
The abundances of Behr (2003) were derived from data with 
a resolution of 0.1 A, while the data used in this paper have 
a resolution of about 0.7 A. Between 11000K and 13 500 K 
the abundances show a clear correlation with the resolution of 
the data from which they were derived - the abundances de- 
crease with increasing resolving power. This behaviour had 
been noticed already by Moni Bidin et al. (2012. This may be 
caused by the fact that at these relatively cool temperatures the 
helium lines are rather weak, which together with the helium 
deficiency makes them hard to fit. For higher temperatures the 
abundances from the different sources overlap, with the excep- 


15000 10000 8000 



Fig. 6. Masses from line profile fits. For comparison we also show a 
canonical zero-age horizontal branch from Moehler et al. (2003). In 
the upper plot we show also results from FORS2 observations of hot 
HB stars in M 80, NGC 5986 (Moni Bidin et al. 2009), NGC 6752 
(Moni Bidin et al. 2007), and M22 (Salgado et al. 2013). In the lower 
plot we also provide the results obtained from FLAMES and FORS2 
observations of hot horizontal branch stars in wCen (Moehler et al. 
2011, small squares; Moni Bidin et al. 2011, small triangles). 


tion of the hottest star of Behr (2003). Due to the small number 
of results from high-resolution spectroscopy it is not clear if the 
resolution effect persists to higher temperatures. For a compar- 
ison to the predictions by diffusion theory see Sect. 6.4. 

6.2. Heavy elements 

Although the spectra have only medium resolution we esti- 
mated abundances via spectrum synthesis using the abundance 
fitting routine of LINFOR in order to verify whether our as- 
sumptions about the overall increase in heavy elements was 
roughly correct. We used the line lists from Kurucz (1993) and 
determined simultaneously abundances for He i, Mg n, Si n, P n, 
Tin, Mnn, Fen, and Nin. Once the fitted abundances did not 
change anymore from one iteration to the next, we performed 
a visual check to verify that the spectrum had indeed been 
well reproduced. The results are shown in Fig. 8. The LINFOR 
abundance fitting routine provides only an overall error, but 
the scatter of abundances below 1 1 000 K gives an idea of the 
minimum uncertainties, as we do not expect the abundances in 
shown in Fig. 8 to vary in the cool stars. For the stars hotter than 
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Fig. 7. Helium abundances for hot HB stars in M80, NGC 5986, 
NGC 6752, and tuCen (Moni Bidin et al.), M 3, M 13, and NGC 6752 
(Moehler et al.) and NGC 288. The grey dots mark the abundances 
predicted for various HB ages from diffusion theory (with an ad hoc 
surface mixing zone, Michaud et al. 2011, see Sect. 6 for details). The 
tracks have dots every 3 Myr. 


the u - jump additional uncertainties come from the fact that we 
use homogeneous model atmospheres to analyse stars affected 
by diffusion. 

The results for P n, Ti n, Mn n, Fe n, and Ni n agree roughly 
with our assumption of [M/H] = +0.5 for the stars hotter than 
the w-jump, while Mgn and Sin, on the other hand, are rather 
lower than our assumed abundances. However, since the atmo- 
spheric structure depends more on iron than on the light ele- 
ments we take these results as indication that our assumptions 
were reasonable. For the cool stars, Mg n and Ti n seem to show 
an unexpected gradient with temperature, but it is unclear if this 
is significant. 

Our results show some differences compared to those ob- 
tained for NGC 288 by Behr (2003). Behr (2003) found no clear 
dependence of the Mg abundance relative to T e ff for the three 
clusters studied with a sufficient number of stars with T e g above 
the w-jump, while a trend appears to be present in our results for 
NGC 288. For Si and Mn Behr found significantly lower abun- 
dances than we do, which might be an effect of instrumental 
resolution (see Sect. 6.1). 

6.3. Abundances predicted by stratified model spectra 

In order to compare the results from our observed spectra to 
the predictions from stratification theory we also determined 
abundances from model spectra of LeBlanc et al. (2009) that 
include abundance stratification due to diffusion. These mod- 
els calculate self-consistently the structure of the atmosphere 
with the vertical abundance stratification. The abundances of 
the individual elements are calculated in each layer of the at- 
mosphere while assuming equilibrium (a nil diffusion veloc- 
ity). This leads to vertical abundance stratification and modi- 
fies the atmospheric structure. First we compared the observed 
metal lines to those predicted by diffusion models and found 
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Fig. 8. Abundances derived via spectrum synthesis for all stars hotter 
than 9 000K. The triangles mark the same stars as in Fig. 1. The as- 
terisks mark the results from Behr (2003). The bars at 17 500 K mark 
the average error bars of Behr (2003) for NGC 288. Our errors are 
unlikely to be smaller. The dashed lines mark the solar abundances. 
The 4-pointed stars mark the abundances derived for stratified model 
spectra (left column, see Sect. 6 for details), which give an indication 
of the equilibrium abundances achievable by diffusion. The grey dots 
(right column) mark the abundances predicted for various HB ages 
from diffusion theory (with an ad hoc surface mixing zone, Michaud 
et al. 2011, see Sect. 6.4 for details). The tracks have dots every 1 Myr. 


that the observed lines were much weaker than predicted by 
these models (see Fig. 9). 

As discussed in LeBlanc et al. (2009), the abundances of 
some elements (Fe for example) predicted by the models can 
be overestimated because in a real star the diffusion process is a 
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Table 5. Temperatures and surface gravities obtained from fitting strat- 
ified model spectra with homogeneous model spectra for [M/H] =+l 
(see Sect. 6 for details). The errors are derived from the )f fit (see 
Sect. 5.2 for details) and are statistical errors only. 



Model 

Teff 

logg 

Fit 

Teff 

logg 

A 

12000 

3.5 

1 1 500± 1 00 

3.07±0.05 

B 

12000 

4.0 

1 1 700± 1 00 

3.60±0.06 

C 

12000 

4.5 

1 1 300± 1 00 

3.94±0.06 

D 

13000 

3.5 

12200±100 

3.08±0.05 

E 

13000 

4.0 

12000±100 

3.48±0.05 

F 

13000 

4.5 

12400±100 

4.00±0.05 

G 

14000 

4.0 

12700±100 

3.47±0.04 

H 

14000 

4.5 

13200±100 

4.01 ±0.04 

I 

16000 

4.5 

14300±100 

3.98±0.04 

J 

18000 

5.0 

16300±100 

4.52±0.03 


time-dependent phenomenon. Even though the radiative forces 
on a given element can theoretically support a very large over- 
abundance, this situation may be hampered by various factors. 
The equilibrium solution used in these models may therefore 
not be reached. For instance, in a real star, such a large quantity 
of atoms might not be able to come to the surface due to the 
diffusion that takes place below the stellar atmosphere. Also, 
the abundances of other elements (such as helium) which have 
relatively weak radiative accelerations can be underestimated 
by these models. 

Next we convolved the stratified model spectra to the same 
resolution as our observed data and multiplied them with a 
spectrum of average 1 and rms of 0.0125, resulting in a signal- 
to-noise ratio of 80. Then we determined their effective tem- 
peratures and surface gravities in the same way as described 
in Sect. 5.2, except that we used a model grid with [M/H] = 
+1, since the lines in the stratified model spectra were much 
stronger than in the observed spectra (see Fig. 9). The resulting 
parameters are given in Table 5 and clearly show that the analy- 
sis of stratified model spectra with homogeneous model spectra 
yields lower temperatures and surface gravities, with the differ- 
ence in temperature increasing with increasing temperature. 

Using the parameters in the right column of Table 5 we 
then determined abundances for those model spectra with fit- 
ted log g above 3.45 as expected for stars on or near the ZAHB 
at the fitted temperatures, i.e. all models except A and D. We 
alternated between simultaneously fitting the ions Mg n, Si n, 
P n, S n, A n, Ti n, V n, Cr n, and simultaneously fitting the ions 
Mnn, Fei, Fen, Con, Nin, Srn, Zrn, until the fit could not 
be improved. Then we checked again visually that most of the 
spectrum had been reproduced. Usually we found that some 
lines were not fit well, but that is to be expected since the 
abundance distribution in the model atmospheres creating these 
spectra is far from homogeneous, resulting in a very different 
atmospheric structure. The results of these fits are marked by 
four-pointed stars in the left column of Fig. 8. As already sug- 
gested by the comparison shown in Fig. 9 the abundances de- 
rived from the stratified spectra are generally higher than those 
derived from the observed spectra, with He i (virtually absent 


from stratified model atmospheres) and Mgn being the excep- 
tions. 

6.4. Abundances predicted by diffusion models 

It is also possible to compare our abundance results to stellar 
evolution models computed including the effects of atomic dif- 
fusion (Michaud et al 2008, 2011). These models followed the 
evolution from the zero-age main sequence treating in detail 
atomic diffusion inside the star. On the horizontal branch, they 
predicted large overabundances of metals, often larger than ob- 
served above 1 1 000 K (similar to the problems found above 
for the stratified model atmospheres). A surface mixed mass of 
around 1 0 7 M Q was needed to reduce the expected anomalies 
to values observed in HB stars of a number of clusters as well 
as in sdB stars. However, the introduction of a surface mix- 
ing layer leads to a vertically homogeneous atmosphere. This 
is in contradiction to the observed vertical stratification of cer- 
tain elements, including iron, in some blue HB stars (Khalack 
et al. 2007, 2008, 2010). The tracks shown in Fig. 8 are from 
Michaud et al. (2011). Those tracks start with the original 
cluster abundances at the ZAHB and cover the first third of 
horizontal-branch evolution which lasts some 100 Myr 2 . 

Fig. 8 shows that the abundances change more rapidly dur- 
ing the first 10 Myr than during the following 20 Myr for all 
species that become overabundant. This is largely caused by 
a reduction of the radiative accelerations as the concentrations 
increase. The abundance increase is also expected to be slow 
during the following 70 Myr for those species. The Si, Ti, Fe, 
and Ni observations are within the expected abundance range 
when one takes error bars evaluated from the scatter of Fe 
and Si below 1 1 000 K into account. For Si the trend sug- 
gested by the models appears to be present in the data. For 
Mn our abundances are higher than the results from Behr et 
al. (2003), which agree well with the predictions. Phosphorus 
is some 5 times more overabundant than predicted by the mod- 
els, whereas magnesium shows small effects while the models 
predict larger abundances by a factor of order 5. 

Helium is observed to be more underabundant than ex- 
pected after 30 Myr (cf. Fig. 7). However its abundance has a 
very different time dependence from that of metals. The under- 
abundance of helium is caused by gravitational settling and this 
process is as rapid from 10 to 30 Myr as during the first 1 0 Myr. 
Since helium settles, the settling goes as e - '^ where 9 is the set- 
tling time scale. So the settling continues exponentially and if 
helium is underabundant by a factor of 3 after 30 Myr it should 
be underabundant by a factor of around 27 after 100 Myr, corre- 
sponding to a value of 9.56 in Fig. 8 at the end of the horizontal 
branch evolution. Since it is highly unlikely, however, that most 
of the stars in Fig. 8 are at the end of their HB evolution, some 
discrepancy remains. 

As may be seen from Fig. 4 of Michaud et al (201 1), abun- 
dance anomalies caused by atomic diffusion are not limited to 
a thin surface phenomenon. However, there are species (e.g. 

2 The models around 12 000 and 13 500K were stopped around 
10 Myr because of convergence problems, but the trend can be esti- 
mated from the surrounding models. 
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Mg) whose observed anomalies do not agree with the expected 
ones suggesting that the model may be missing something. 
Assuming that the outer 1 0 7 M Q is completely mixed could 
be an oversimplification, especially since the mixing mecha- 
nism is currently unknown. It could for instance be thermoha- 
line convection (see Sect 5.3 of Michaud et al 2011). This is rel- 
atively weak convection which might not eliminate all effects 
of additional diffusion in the atmosphere. For instance, helium 
is largely neutral in the atmosphere of HB stars of 1 1 000 to 
15 000 K so that it has an atomic diffusion coefficient larger 
than that of ionized metals (by a factor of around 100); helium 
might then be the most affected atomic species as turbulence 
weakens. This remains speculative. 



4500 4520 4540 4560 4580 4600 

X[A] 

Fig. 9. Nonnalised spectra for stars 199 (upper panel, 12400K) and 
221 (lower panel, 13 400K) compared to stratified model spectra F 
(upper panel) and H (lower panel), which have fitted parameters close 
to the ones of the observed spectra. For clarity, the stratified model 
spectra have been offset by 0.1 from the observed spectra (see Tables 4 
and 5). 


7. Line profile fitting with stratified model spectra 

We also verified how the use of stratified instead of homoge- 
neous model atmospheres to fit our observations influences the 
parameters determined from these fits. In order to allow a di- 
rect comparison to results obtained using homogeneous model 
spectra without metal lines to fit our observations, we calcu- 
lated a set of synthetic spectra with only H present (the abun- 
dance of all of the other elements were set at -5.00 while H is 


at 12.00) while using the relations between temperature, den- 
sity, and radius from the stratified model atmospheres. Helium 
lines were not included as the stratification abundances of he- 
lium are extremely low. For comparison we also fitted the ob- 
served spectra again (H lines only) with homogeneous model 
spectra computed with [M/H] = +0.5 and log ^ = -3. The 
results are shown in Fig. 10. The prediction of larger spectro- 
scopic gravities by the models including abundance stratifica- 
tion shown here was also found in Hui-Bon-Hoa et al. (2000) 
and LeBlanc et al. (2010). The amplitude of the effect might be 
overestimated by the stratified models due to the fact that, in 
most cases, they predict larger abundance anomalies for blue 
HB stars than observed (see Sect. 6). 



Fig. 10. Atmospheric parameters derived from the line profile fitting 
for the stars hotter than the (/-jump, using homogeneous model at- 
mospheres ([M/H] = +0.5, log ^ = -3, filled circles) and stratified 
model atmospheres (filled triangles). For comparison we also show 
a canonical zero-age horizontal branch and temiinal age horizontal 
branch from Moehler et al. (2003). 

In Fig. 10 all stars move to higher temperatures and sur- 
face gravities when fitted with model spectra from stratified 
model atmospheres (as expected from the results of the reverse 
fitting given in Table 5), but only the stars between 14 000 K 
and 16 000K move actually closer to the ZAHB. Cooler stars 
show too high surface gravities, while hotter ones move paral- 
lel to the TAHB. The principal defect of the stratified models is 
that they generally predict larger abundances than observed. In 
more realistic models, one would expect that the shifts shown 
in Fig. 10 be somewhat smaller. The lower temperature stars 
could then possibly fall between the TAHB and ZAHB in the 
logg - Ten plane thus giving satisfactory results. 

As a consistency check we derived masses for the new 
parameters, adjusting the theoretical y-magnitudes of Kurucz 
(1993) for [M/H] = +0.5 by -0T3 to account for the strati- 
fication effects. This offset was determined by comparing the 
fluxes of stratified and homogeneous model spectra. The results 
are shown in Fig. 11. Here the values obtained from stratified 
model spectra for stars above 14 000 K are now closer to the 
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canonical values than those obtained from homogeneous model 
spectra. This provides some support to the notion that the too 
low surface gravities and masses found from analyses with ho- 
mogeneous model spectra are due to the mismatch in atmo- 
spheric structure between homogeneous and stratified model 
atmospheres. One should keep in mind, however, that the strat- 
ified model atmospheres predict in general much stronger lines 
than are observed, so it is currently not clear if a fully self- 
consistent solution can be achieved. Similarly to the discussion 
related to the large shifts found in Fig. 10, the masses predicted 
in Fig. 11 by more realistic stratified models for stars just above 
the temperature threshold of the w-jump, in which the abun- 
dances are less extreme, would be closer to the ZAHB. 



Fig. 11. The masses detemiined from line profile fits for the stars bluer 
than the Grundahl jump, using homogeneous (filled circles) and strat- 
ified (filled triangles) model spectra. 


8. Overluminous Stars 

One possible explanation for the overluminous stars mentioned 
in Sect. 1.2 is that they have evolved from hotter locations blue- 
ward of the w-jump and are now near the end of their HB phase 
evolving towards the asymptotic giant branch. As the evolu- 
tion away from the zero-age HB increases the overall luminos- 
ity of an HB star, one would expect such evolved stars to be 
overluminous in other bandpasses besides u compared to stars 
near the zero-age HB. Fig. 12 shows the relation between the 
effective temperature derived from line profile fits and the y 
magnitude of the HB stars. Here especially the hot overlumi- 
nous stars are clearly brighter than the ’normal’ stars at similar 
temperatures. To a lesser extent this is also true for the cool 
overluminous stars. In Figs. 4 and 5, however, only three of the 
overluminous stars, namely 100 (hot), 103, and 127 (both cool) 
are clearly separated from the majority of the stars by a signif- 
icantly lower surface gravity. It is not obvious why the other 
overluminous stars have higher luminosities, as they are incon- 
spicuous with respect to the majority of the stars in all other 
parameters (effective temperature, surface gravity, and abun- 
dances). However, even three stars evolving off the ZAHB in 


the temperature range 9 700 K to 11 500 K pose a problem for 
evolutionary time scales, as one would expect to find of the or- 
der of 100 HB stars close to the ZAHB between about 1 1 000 K 
and 14 000 K (corresponding to a range of roughly 0.5 - 1.05 in 
u —y in Fig. 1) for each of the evolved stars, which are clearly 
not seen. 
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Fig. 12. The y magnitudes of the HB stars compared to the effective 
temperatures derived from line profile fits. 

If the overluminous stars are post-HB stars evolved from a 
hotter location along the blue tail (as suggested by their higher 
luminosity), one would expect them to preserve their abun- 
dance anomalies as they approach the w-jump until dredge-up 
from the deepening hydrogen convection zone causes the sur- 
face abundances to return approximately to the normal cluster 
abundances. However such stars are expected to have had a 
low He abundance on the HB and a low rotation rate which 
may cause them to maintain their abundance anomalies to tem- 
peratures somewhat lower than the w-jump. First, due to the 
depletion of their surface helium, the overluminous stars just 
redward of the w-jump may have slightly shallower convection 
zones than the other stars of the same effective temperature. 
However hydrogen ionization is the main driver of the con- 
vection when the effective temperature is below 10 000 K and 
it will lead to progressive helium dredge-up so that the differ- 
ence in convection zone depth could be small below 10 000 K. 
Second, due to the rotation rate drop in HB stars blueward of 
the u-jump, the overluminous stars would be expected to have 
low rotation rates, if they are evolved from the blue HB. This, 
in turn, would reduce the meridional flows which inhibit diffu- 
sion in stars cooler than the w-jump. However, as may be seen 
from Fig. 4 of Quievy et al. (2009), the limiting rotation rate 
decreases very rapidly as the effective temperature drops below 
10 000 K; a star with a v e = 8 km s 1 would show large effects of 
diffusion above 12 000 K but would be mixed below 10 000 K. 
The results given in Fig. 8 indeed show that the effects of radia- 
tive diffusion disappear between the overluminous stars hotter 
than the w-jump at 1 1 500 K and the overluminous stars cooler 
than the w-jump at 10 000 K. So it would appear that the low 
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He abundance due to diffusion and low rotation do not prevent 
the effects of diffusion from disappearing rather quickly as a 
star evolves through the zz-jump to cooler temperatures. This 
result supports the general assumption that the abundances of 
HB stars redward of the zz-jump are not affected by diffusion, 
irrespective of their evolutionary status. The observations are 
ahead of the stellar evolution models in this case. To put the 
theory on a stronger base, one would need to continue the 
evolution calculations with diffusion to the end of the HB; 
however one would also need to include meridional circula- 
tion to verify at what effective temperature the mixing oc- 
curs and how this temperature depends on rotation. 


9. Summary and Conclusions 

- The atmospheric parameters and masses of the hot HB stars 
in NGC 288 show a behaviour seen also in other clusters for 
temperatures between 9 000K and 14 000 K. Outside this 
temperature range, however, they follow the results found 
for such stars in a> Cen. 

- The abundances derived from our observations are for most 
elements within the abundance range expected from evo- 
lutionary models that include the effects of atomic diffu- 
sion and assume a surface mixed mass of 10 7 M Q as deter- 
mined previously for sdB stars and other clusters. The ex- 
ceptions are helium, which is more deficient than expected 
and phosphorus, which is substantially more abundant than 
predicted. The abundances predicted by stratified model at- 
mospheres, which were not adjusted to observations, are 
generally significantly more extreme than observed, except 
for magnesium. 

- Analysing the observed spectra with stratified model spec- 
tra moves the HB stars to higher temperatures, surface grav- 
ities, and masses. Since the equilibrium abundances lead to 
excessive adjustments, more realistic abundance gradients 
may well lead to models that put the HB stars between the 
TAHB and ZAHB (see Sect. 7). Model atmospheres includ- 
ing such improvements are needed to answer this question. 

- Five of the eight overluminous stars around the zz-jump, 
which we observed, do not deviate substantially from the 
other HB stars in the same temperature range in any of the 
parameters we determined. We are therefore at a loss to ex- 
plain their brighter luminosities. The remaining three over- 
luminous stars do show lower surface gravities, as would be 
expected if they evolve away from the HB. However, they 
pose a substantial problem for evolutionary time scales, as 
one would expect to find of the order of 100 HB stars close 
to the ZAHB between about 1 1 000 K and 14 000 K (corre- 
sponding to a range of roughly 0.5 - 1.05 in u -y in Fig. 1) 
for each of the evolved stars, which are clearly not seen. 

- With respect to their abundances all overluminous stars 
show the same behaviour as the majority of the stars in 
the respective temperature range. This result supports the 
general assumption that the abundances of HB stars 
redward of the z/-jump are not affected by diffusion, ir- 
respective of their evolutionary status. 


The evolution models including diffusion and the stratified 
model atmospheres both predict higher than observed abun- 
dances for many elements affected by radiative levitation. The 
evolution models can be adjusted to reproduce the observed 
abundances by introducing an ad hoc defined mixed zone, 
which, however, is potentially inconsistent with the observed 
vertical stratification of at least some elements in the atmo- 
sphere. For the stratified model atmospheres, it is currently 
unclear whether a more limited abundance stratification which 
would provide a better description of the observed abundances 
could still explain the photometric anomalies around the zz- 
jump. Our results show definite promise towards solving the 
long-standing issue of surface gravity and mass discrepancies 
for hot HB stars, but there is still much work needed to arrive 
at a self-consistent solution. 
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Table 1 . Target coordinates and photometric data. 


ID 

® 2000 

<52000 

y 

b-y 

Cl 

m. 

52 

00:52:50.77 

-26:38:02.9 

15.357±0.001 

+0.030±0.003 

+ 1.301±0.005 

+0.131±0.005 

55 

00:52:47.43 

-26:33:14.2 

15.318±0.002 

+0.094±0.002 

+ 1.253±0.004 

+0.122±0.004 

60 

00:52:38.59 

-26:37:05.1 

15.328±0.003 

+0.110±0.004 

+ 1.202±0.008 

+0.118±0.007 

61 

00:52:42.80 

-26:36:39.2 

15.3 12±0.001 

+0.127±0.002 

+ 1.151 ±0.004 

+0.121±0.003 

63 

00:52:47.08 

-26:35:25.0 

15.345±0.002 

+0.094 ±0.003 

+ 1.264±0.005 

+0.119±0.005 

70 

00:52:45.40 

-26:35:21.5 

15.438±0.002 

+0.055±0.003 

+ 1.288±0.004 

+0.133±0.004 

72 

00:52:53.92 

-26:38:45.6 

15.421±0.002 

+0.088±0.004 

+ 1.234±0.008 

+0.125±0.008 

74 

00:52:35.47 

-26:34:24.8 

15.453±0.001 

+0.055±0.002 

+ 1.283±0.005 

+0.133±0.003 

79 

00:52:42.68 

-26:34:50.2 

15.507±0.002 

+0.050±0.003 

+ 1.307±0.004 

+0.135±0.005 

81 

00:52:40.77 

-26:33:47.7 

15.569±0.001 

+0.034±0.002 

+ 1.288±0.003 

+0.137±0.003 

83 

00:52:32.07 

-26:35:46.6 

15.560±0.002 

+0.022±0.002 

+ 1.266±0.007 

+0.139±0.004 

86 

00:52:52.51 

-26:34:29.7 

15.582±0.001 

+0.038±0.002 

+ 1.284±0.006 

+0.136±0.004 

88 

00:52:52.10 

-26:34:12.1 

15.639±0.002 

+0.021 ±0.003 

+ 1.260±0.007 

+0.131±0.005 

89 

00:52:37.88 

-26:36:35.4 

15.518±0.002 

+0.042±0.002 

+ 1.242±0.006 

+0.134±0.003 

90 

00:52:46.64 

-26:39:03.1 

15.634±0.002 

+0.020±0.004 

+ 1.253±0.011 

+0.131±0.009 

96 

00:52:37.64 

-26:31:13.5 

15.735±0.004 

+0.005±0.004 

+ 1.238±0.010 

+0.134±0.007 

99 

00:52:29.92 

-26:36:07.4 

15.799±0.002 

-0.004±0.003 

+ 1.182±0.009 

+0.133±0.007 

100 

00:52:38.68 

-26:35:58.6 

15.959±0.002 

-0.034±0.002 

+0.858±0.004 

+0.112±0.004 

101 

00:53:04.16 

-26:38:29.8 

15.91 1±0.005 

-0.024±0.009 

+0.938±0.01 1 

+0.11 6±0.0 1 5 

102 

00:52:51.47 

-26:36:26.0 

1 5 .8 1 5±0.001 

+0.007±0.002 

+ 1.205±0.004 

+0.128±0.003 

103 

00:52:35.45 

-26:39:11.0 

15.840±0.003 

-0.01 1±0.004 

+ 1.063±0.010 

+0.115±0.007 

106 

00:53:02.55 

-26:35:32.9 

15.822±0.002 

+0.005±0.003 

+ 1.176±0.008 

+0.127±0.005 

107 

00:52:32.39 

-26:36:30.2 

15.810±0.002 

+0.008±0.002 

+ 1.173±0.009 

+0.126±0.004 

111 

00:53:05.31 

-26:32:45.2 

15.982±0.003 

-0.017±0.004 

+ 1.084±0.005 

+0.136±0.005 

113 

00:52:39.48 

-26:36:45.7 

15.866±0.001 

+0.001 ±0.002 

+ 1.161 ±0.006 

+0.125±0.004 

114 

00:52:38.68 

-26:37:49.6 

15.905±0.004 

-0.002±0.006 

+ 1.132±0.015 

+0. 126±0.0 1 1 

115 

00:52:45.24 

-26:37:55.1 

1 5.88 1 ±0.001 

+0.000±0.002 

+ 1.170±0.006 

+0.125±0.005 

118 

00:52:55.50 

-26:35:08.2 

15.902±0.001 

-0.005±0.002 

+ 1.158±0.008 

+0.124±0.003 

119 

00:52:56.84 

-26:33:44.8 

15.883±0.001 

+0.002±0.002 

+ 1.159±0.007 

+0.123±0.004 

120 

00:53:01.85 

-26:37:53.2 

15.913±0.004 

-0.008±0.006 

+ 1.136±0.007 

+0.128±0.010 

122 

00:52:37.78 

-26:39:31.6 

15.880±0.002 

-0.012±0.004 

+ 1.094±0.011 

+0.127±0.009 

127 

00:52:39.32 

-26:34:31.6 

15.953±0.002 

-0.014±0.002 

+ 1.047±0.004 

+0.1 1 1±0.004 

142 

00:52:50.55 

-26:36:49.8 

16.061±0.002 

-0.032±0.002 

+0.876±0.005 

+0.1 1 1±0.004 

143 

00:52:52.77 

-26:34:53.0 

16.003±0.002 

-0.005±0.003 

+ 1.084±0.006 

+0.119±0.006 

145 

00:52:40.58 

-26:32:48.6 

15.981±0.002 

-0.013±0.002 

+ 1.108±0.004 

+0.124±0.003 

146 

00:52:43.35 

-26:37:56.4 

16.098±0.002 

-0.019±0.004 

+0.878±0.005 

+0.107±0.006 

147 

00:52:37.26 

-26:36:46.9 

16.033±0.002 

-0.013±0.003 

+ 1.054±0.007 

+0.127±0.004 

149 

00:52:33.14 

-26:33:44.6 

16.084±0.003 

-0.020±0.004 

+ 1.044±0.009 

+0.122±0.005 

151 

00:52:48.32 

-26:32:57.6 

16.032±0.002 

-0.016±0.003 

+ 1.074±0.004 

+0.131±0.004 

154 

00:52:54.89 

-26:37:11.7 

16.058±0.001 

-0.007±0.002 

+ 1.054±0.007 

+0.112±0.004 

156 

00:52:50.53 

-26:35:12.0 

16.039±0.003 

-0.004±0.005 

+ 1.060±0.007 

+0.112±0.009 

157 

00:52:53.76 

-26:39:08.7 

15.992±0.004 

-0.013±0.007 

+ 1.073 ±0.011 

+0.123±0.013 

167 

00:52:46.42 

-26:34:07.7 

16.093±0.002 

-0.009±0.003 

+ 1.039±0.006 

+0.112±0.006 

169 

00:52:46.50 

-26:31:30.7 

16.153±0.002 

-0.023±0.003 

+0.957±0.006 

+0.119±0.005 

176 

00:52:48.70 

-26:34:00.7 

16.140±0.002 

-0.012±0.003 

+0.984±0.006 

+0.113±0.006 

179 

00:52:48.17 

-26:35:19.9 

16.236±0.003 

-0.032±0.006 

+0.844±0.009 

+0.108±0.01 1 

180 

00:52:50.91 

-26:36:09.4 

16.153±0.002 

-0.014±0.003 

+0.975±0.004 

+0. Ill ±0.004 

183 

00:52:39.91 

-26:37:23.8 

16.254±0.002 

-0.025±0.004 

+0.784±0.01 1 

+0.105±0.009 

187 

00:52:44.71 

-26:35:31.4 

16.307±0.003 

-0.033±0.005 

+0.843±0.008 

+0.116±0.009 

195 

00:52:27.41 

-26:35:58.8 

16.422±0.004 

-0.050±0.005 

+0.685±0.008 

+0.122±0.008 

196 

00:52:44.29 

-26:35:53.2 

16.357±0.002 

-0.031 ±0.002 

+0.760±0.004 

+0.105±0.004 

199 

00:52:55.57 

-26:32:58.7 

16.425±0.002 

-0.042±0.002 

+0.726±0.005 

+0.113±0.004 

212 

00:52:59.33 

-26:39:00.4 

16.605±0.005 

-0.068±0.009 

+0.577±0.013 

+0. 135±0.0 1 7 

213 

00:52:42.83 

-26:31:06.8 

16.627±0.005 

-0.062±0.005 

+0.564±0.008 

+0.128±0.007 
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Table 1. Target coordinates and photometric data (cont’d) 


ID 

® 2000 

<52000 

y 

b-y 

Cl 

m. 

216 

00:52:54.57 

-26:33:20.4 

16.572±0.002 

-0.051 ±0.005 

+0.627±0.01 1 

+-0. 1 17 ±0.0 1 0 

221 

00:52:52.31 

-26:35:13.7 

16.638±0.002 

-0.037±0.003 

+0.584±0.007 

+0.105±0.006 

228 

00:52:47.36 

-26:37:52.6 

16.761±0.002 

-0.057±0.004 

+0.495±0.009 

+0.116±0.008 

230 

00:52:24.33 

-26:35:23.5 

16.759±0.009 

-0.059±0.014 

+0.485±0.029 

+0.114±0.025 

231 

00:52:47.01 

-26:36:23.0 

16.714±0.003 

-0.057±0.004 

+0.552±0.008 

+0.115±0.006 

240 

00:52:43.69 

-26:35:01.7 

16.761±0.003 

-0.053±0.005 

+0.534±0.007 

+0.113±0.008 

242 

00:52:45.02 

-26:37:35.4 

16.828±0.002 

-0.064±0.005 

+0.421 ±0.007 

+0.121 ±0.009 

243 

00:52:44.02 

-26:35:42.2 

16.889±0.003 

-0.043±0.003 

+0.445±0.006 

+0.105±0.005 

275 

00:52:49.39 

-26:35:53.7 

17.030±0.006 

-0.045±0.009 

+0.391±0.015 

+0.098±0.016 

288 

00:52:49.30 

-26:38:19.1 

17.1 15±0.003 

-0.058±0.004 

+0.358±0.006 

+0.104±0.006 

292 

00:53:00.03 

-26:36:32.9 

17.1 17±0.003 

-0.068±0.005 

+0.406±0.008 

+0.108±0.008 

300 

00:52:49.11 

-26:35:35.1 

17.106±0.003 

-0.048±0.007 

+0.362±0.009 

+0.100±0.013 

304 

00:52:48.60 

-26:33:17.1 

17.173±0.003 

-0.054±0.005 

+0.307±0.01 1 

+0.107±0.009 

347 

00:52:50.31 

-26:38:30.0 

17.532±0.003 

-0.066±0.005 

+0.275±0.009 

+0.102±0.009 


Table 2. Observing times, conditions and setups 


Date 

Start 

[UT] 

Seeing 

Airmass 

Moon 

Ilium. Dist. 

setup 

exposure 
time [s] 

2003-07-04 

08:15:23 

175 

1.15 

0.21 

14619 

LR1 

2520 

2003-07-04 

09:06:06 

171 

1.06 

0.22 

14713 

LR1 

2520 

2003-07-04 

09:49:14 

171 

1.02 

0.22 

14717 

LR1 

2000 

2003-07-05* 

07:11:48 

170 

1.37 

0.30 

15312 

LR1 

2520 

2003-07-05* 

07:56:47 

172 

1.21 

0.30 

15314 

LR1 

1432 

2003-07-07 

08:04:26 

079 

1.15 

0.52 

14917 

LR1 

2520 

2003-07-07 

08:49:35 

078 

1.06 

0.53 

14914 

LR1 

2520 

2003-07-07 

09:33:38 

078 

1.02 

0.53 

14912 

LR1 

2520 

2003-07-07* 

10:16:53 

078 

1.00 

0.53 

14819 

LR1 

832 

2003-07-08 

08:23:19 

076 

1.10 

0.63 

14012 

LR1 

2520 

2003-07-08 

09:09:37 

077 

1.03 

0.64 

13918 

LR1 

2520 

2003-07-08* 

09:53:34 

078 

1.01 

0.64 

13914 

LR1 

1417 

2003-07-09 

09:19:31 

0 77 

1.02 

0.75 

12812 

LR1 

2520 

2003-07-10 

09:06:35 

076 

1.03 

0.84 

11519 

LR3 

3000 

2003-07-10 

10:01:51 

171 

1.00 

0.84 

11513 

LR3 

661 

2003-07-20 

05:36:53 

077 

1.59 

0.59 

3011 

LR3 

2520 

2003-07-20 

06:20:22 

170 

1.40 

0.59 

3013 

LR3 

623 

2003-07-21* 

06:15:02 

076 

1.34 

0.50 

3810 

LR3 

2520 

2003-07-21 

06:58:04 

078 

1.18 

0.49 

3812 

LR3 

2520 

2003-07-21 

07:47:32 

172 

1.07 

0.49 

3815 

LR2 

2520 

2003-07-21 

08:30:33 

175 

1.02 

0.49 

3817 

LR2 

2520 

2003-07-21 

10:17:49 

171 

1.02 

0.48 

3913 

LR2 

1736 

2003-07-22 

05:03:19 

074 

1.83 

0.41 

4616 

LR1 

2520 

2003-07-22 

05:46:25 

074 

1.47 

0.40 

4619 

LR1 

2520 

2003-07-22 

06:37:10 

075 

1.24 

0.40 

4713 

LR2 

2520 

2003-07-22 

07:22:53 

075 

1.11 

0.40 

4716 

LR2 

2520 

2003-07-22 

08:12:00 

075 

1.04 

0.40 

4719 

LR2 

2520 

2003-07-22* 

10:21:07 

076 

1.02 

0.39 

4816 

LR2 

20 

2003-07-27 

07:45:34 

079 

1.04 

0.04 

10112 

LR1 

2520 

2003-07-28 

05:59:11 

075 

1.29 

0.01 

11113 

LR2 

2520 

2003-07-28 

06:49:13 

075 

1.13 

0.01 

11119 

LR1 

2520 


*: not used in final analysis, usually due to S/N 
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Table 3. Temperatures and surface gravities for stars redder than the j(-jump and not showing any evidence of radiative levitation as derived 
from the a 0 and sz parameters and from line profile fits for an assumed metallicity of [M/H] = —1.5. The errors for the photometric results are 
derived from the errors of a o and sz in combination with the errors of the fits and an uncertainty of 0™02 for the reddening (see Sect. 5.1 for 
details). The errors for the spectroscopic results are derived from the^ 2 fit (see Sect. 5.2 for details) and are statistical errors only. In addition 
we list the helium abundance as determined with LINFOR (see Sect. 6 for details). 


ID 

ao 

[mmag] 

Photometric Parameters 
sz T t ff 

[mmag] [K] 

logg 

Spectroscopic Parameters 
^eff logg log(^) L INFOR 

[K] 

52 

+48± 

4 

91± 5 

9100+380 

3.39+0.11 

8240+ 20 

2.87+0.02 


55 

+ 124± 

4 

91± 4 

8100+260 

3.08+0.15 




60 

+ 134± 

6 

89± 7 

8000+260 

3.05+0.19 




61 

+ 151+ 

3 

79± 3 

7900+210 

3.09+0.16 




63 

+ 125± 

5 

95± 5 

8100+270 

3.02+0.17 




70 

+82± 

4 

86± 4 

8600+330 

3.36+0.12 

7990+ 20 

2.94+0.03 


72 

+ 1 14± 

7 

86± 8 

8200+300 

3.21+0.18 




74 

+81+ 

3 

85± 3 

8600+320 

3.37+0.11 

7910+ 10 

2.94+0.04 


79 

+79± 

5 

86± 5 

8600+340 

3.37+0.13 

7990+ 10 

2.97+0.03 


81 

+54± 

3 

84± 3 

9000+350 

3.47+0.10 

8210+ 20 

2.98+0.02 


83 

+34± 

4 

80± 4 

9300+400 

3.57+0.10 

8360+ 30 

3.02+0.02 


86 

+57± 

3 

84 ± 4 

8900+360 

3.46+0.10 

8160+ 20 

3.00+0.02 


88 

+30± 

5 

88± 6 

9300+420 

3.46+0.11 

8390+ 30 

2.98+0.02 


89 

+56± 

3 

81± 3 

8900+360 

3.51+0.10 

8070+ 20 

2.93+0.03 


90 

+27± 

7 

87±10 

9400+470 

3.48+0.14 

8450+ 30 

3.04+0.02 


96 

+5± 

7 

84± 7 

9700+500 

3.54+0.12 

8560+ 30 

3.07+0.02 


99 

-17± 

6 

80± 7 

10100+510 

3.61+0.12 

9130+260 

3.25+0.15 

-0.6 

102 

0± 

3 

86± 3 

9800+440 

3.51+0.09 

9290+320 

3.38+0.18 

-1.4 

103 

-56± 

6 

86± 7 

10900+590 

3.42+0.19 

9850+ 70 

3.34+0.04 

-1.0 

106 

-8± 

5 

84 ± 6 

10000+490 

3.54+0.11 

9230+290 

3.33+0.17 

-1.3 

107 

-5+ 

4 

85± 4 

9900+460 

3.52+0.10 

9350+260 

3.39+0.15 

-1.2 

111 

-52± 

6 

68± 6 

10800+570 

3.82+0.12 

9760+ 90 

3.44+0.05 

-1.0 

113 

-17± 

3 

85± 4 

10100+470 

3.52+0.11 

9210+180 

3.31+0.11 

-0.8 

114 

-27± 

9 

8 1 ± 1 1 

10300+600 

3.58+0.17 

9340+140 

3.33+0.08 

-1.0 

115 

-16± 

4 

86± 5 

10100+480 

3.50+0.12 

9280+180 

3.34+0.10 

-1.1 

118 

-27± 

3 

86± 3 

10300+480 

3.49+0.11 

9440+ 90 

3.38+0.05 

-1.1 

119 

-16± 

3 

87± 4 

10100+470 

3.48+0.11 

9430+170 

3.39+0.09 

-1.2 

120 

-32± 

9 

80±10 

10400+600 

3.60+0.17 

9470+ 70 

3.39+0.04 

-1.0 

122 

-47 ± 

7 

77± 9 

10700+590 

3.64+0.17 

9740+ 80 

3.40+0.04 

-1.1 

127 

-63± 

4 

89± 4 

11000+550 

3.33+0.18 

9920+ 80 

3.37+0.04 

-0.8 

143 

-42± 

5 

84 ± 6 

10600+550 

3.50+0.16 

9810+ 60 

3.45+0.03 

-1.1 

145 

-46± 

3 

82± 3 

10700+520 

3.54+0.12 

9660+ 60 

3.42+0.03 

-1.0 

147 

-54± 

4 

73± 5 

10800+550 

3.71+0.12 

9810+ 80 

3.43+0.04 

-0.8 

149 

-67± 

5 

78± 5 

11100+590 

3.58+0.17 

9800+ 80 

3.41+0.04 

-1.0 

151 

-53± 

4 

72± 4 

10800+540 

3.74+0.11 

9760+ 80 

3.43+0.04 

-0.9 

154 

-51± 

3 

88± 4 

10800+520 

3.39+0.15 

10100+100 

3.51+0.04 

-1.1 

156 

-46± 

8 

88±10 

10700+610 

3.40+0.22 

9760+ 80 

3.37+0.04 

-1.0 

157 

— 52± 1 1 

79±13 

10800+670 

3.59+0.23 

9580+ 90 

3.34+0.05 

-1.1 

169 

-88± 

4 

72± 5 

11600+600 

3.69+0.17 

10400+100 

3.53+0.05 

-1.0 

176 

-72± 

5 

80± 6 

11200+600 

3.52+0.19 

10300+100 

3.52+0.05 

-1.0 

180 

-75+ 

4 

82± 5 

11300+580 

3.46+0.18 

10600+ 80 

3.60+0.04 

-1.1 


